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ABSTRACT 


ill 


General  properties,  including  Hermiticity,  of 
nonlocal  potentials  in  scattering  problems  are  discussed. 

An  alternative  expression  for  flux  density  which  simplifies 
interpretation  of  the  equation  of  continuity  is  proposed. 

Two  attempts  at  fitting  proton-proton  ’S-phase  shift  data 
using  potentials  which  have  nonlocal  character  are  made. 
Examples  for  some  "equivalent"  potentials  are  given.  The 
Klein-Gordon  equation  is  extended  to  include  nonlocal  inter¬ 
actions,  and  the  Green  function  and  Born  approximation  method 
are  used  to  derive  some  relations  of  interest.  Finally, 
the  technique  of  second  quantization  is  applied  to  Schr<5- 
dinger  f  -  field  with  nonlocal  potentials;  it  is  shown 
that  a  many-parti cle  picture  again  emerges,  while  preserving 
the  structure  of  the  Schrddinger  equation. 
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CHAPTER  I.  INTRODUCTION 


The  common  type  of  potentials  used  in  Quantum 
Mechanics  is  one  consisting  of  those  whose  matrices 
<  r  |  V  |  r*>  in  the  r-representation  are  diagonal.  They  are 
known  as  "local"  ( L)  potentials  and  form  a  subset  of  the  set 
of  "nonlocal"  (NL)  potentials,  which  are  those  that  cannot, 
in  general,  be  diagonalized  in  configuration  space. 

We  do  not  have  much  physical  intuitions  about  NL 
potentials.  Moreover,  they  are  mathematically  complicated, 

I! 

Schrodinger  equation  having  in  this  case  the  form  of  an  integro- 
differential  equation.  One  frequently  has  to  resort  to  numeri¬ 
cal  computers.  For  these  reasons,  physicists  have  hesitated 
to  make  use  of  them.  In  the  phenomenological  study  of  nuclear 
forces,  for  example,  L  potentials  have  almost  been  exclusively 
employed . 

Nevertheless,  it  is  becoming  clear  that  some  sorts 
of  nonlocality  must  be  introduced.  There  are  compelling  reasons 
for  doing  so.  As  far  back  as  1936,  Wheeler  ^ ^  established  the 
equivalence  of  the  Majorana  interaction  and  a  certain  simple 
NL  potential.  More  recently,  extensive  studies  of  the  theory 
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( o) 

of  nuclear  matter  by  Bethe  and  Goldstone  and  by  many  others'  ' 
have  shown  that  the  optical  model  potential,  when  derived  as 
the  summation  of  two-body  potentials,  is  nonlocal  in  character. 

The  optical  model  has  been  quite  successful  in  giving 
a  phenomenological  description  of  elastic  scattering  of  nucleons 
by  nuclei.  For  neutrons  of  incident  energy  below  25MeV,  it  has 
been  possible  to  find  a  set  of  optical-model  parameters, 
independent  of  mass  number,  which  vary  smoothly  with  energy(3). 
This  energy  dependence  has  been  explained  by  Perey  and  Buck^^ 
in  terms  of  a  simple  nonlocality  of  the  potential. 

More  specific  considerations  also  lead  one  to  favor 
the  inclusion  of  NL  interactions.  For  example,  as  Breuckner^) 
pointed  out,  the  features  of  the  two-body  interaction  relevant 
to  finite  nuclei  include  (a)  great  strength  of  interaction 
leading  to  nonlocality,  and  (b)  exchange  character  affecting 
nonlocality  of  single-particle  potential.  Recently,  Gammel 
and  Thaler  demonstrated  that  the  proton-proton  scattering 

data  require  that  the  singlet  even  parity  nucleon-nucleon 
interaction  be  nonlocal. 

In  general,  the  potential  must  be  nonlocal  to  represent 
the  effect  of  a  many-body  system  on  the  single  incoming  particle 
(?) .  The  many-body  system  may  be  a  nucleus,  or  the  meson  and 
baryon  cloud  surrounding  a  particle . 
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In  passing,  one  may  mention  that  nonlocality  may 
come  in  in  a  more  fundamental  way,  as  in  Field  Theory . 

*  *  * 

We  give  in  Chap.  II  a  general  discussion  of  NL 
potentials . 

Analytic  solutions  for  the  S-wave  scattering  by 

certain  NL  potential-well  plus  local  parts  are  given  in  Chap.  Ill, 

together  with  attempts  at  fitting  proton-proton  scattering 

phase-shift  data,  for  the  energy  range  20-3^0  MeV. 

In  most  cases,  the  classical  equivalent  of  a  local 

potential  can  be  found  by  inspection.  With  NL  potentials  the 

situation  is  more  involved.  This  and  other  related  topics 

(q) 

will  be  discussed  in  Chap.  IV,  where  the  result  of  Dirac v  ' 
is  used. 

For  high-energy  scatterings,  one  expects  that  the 

M 

Klein-Gordon  equation  would  give  better  results  than  Schrodinger 
equation,  other  things  being  equal.  Some  relevant  relations  for 
NL  interactions  are  given  in  Chap.V. 

In  Chap. VI,  the  technique  of  second  quantization  is 

n 

applied  to  Schrodinger  equation  with  NL  potentials.  The  process 
converts  (as  in  the  local  case)  a  one-particle  theory  to  a  many- 

ii 

particle  theory,  while  preserving  the  structure  of  the  Schrodinger 
equation . 
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CHAPTER  II.  GENERAL  PROPERTIES 


II. 1  Schr8dinger  Equation. 

The  action  of  a  potential  operator  V  on  the  state 

vector  |  f  >  describing  a  single-particle  system  is  represented,, 

—  > 

in  the  r  -  representation.,  by 

<  r  |  V  |  ^  >  —  J  dr  1  <  r  |V|r!>  <  r*\  ip  > 

(II. 1) 

—  J  dr1  V(r,r')  if/  (rT) 

where 

V(r,r’)  -  <r|  V  |  r  !  > 

and 

f  (?)  =  <  r[  if  >  . 

V(r,r ')  is  commonly  referred  to  as  the  "kernel" 

function . 

Correspondingly,  the  time- independent  Schrftdinger 

equation., 

<  r  |  |  ^>  +  <  r  |  V  |  ip  >  =  E  <  r  \  ip  >, 


' 


.  <  ^  I 
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takes  the  form  of  an  integro-dif f erential  equation: 


r 


2m 


V2  +  E)  f( r)  - 


dr  1  V  ( r  „  r. 1 )  ty  ( r  1 ) 


(II. 2) 


\J 


For  local  (L)  potentials.,  i.e.,  when  V(r,r*)  has  the  form 


V(r.,rf)  =  v(r)  5(r-rl) 


(II. 3) 


4- 

0 


he  integration  in  equation  (2)  becomes  trivial.,  and  one  has 


2m 


V2  +  E)  ^(r)  =  v(r)  f{r) . 


II. 2  Angular  and  Energy  Dependence  of  Kernel  Function . Hermiticity . 

Certain  limitations  on  the  kernel  function  V(r,r*)  are 
found  to  be  necessary.  Clearly  they  must  be  such  as  to  ensure 
the  convergence  of  all  quantities  of  physical  interest.  On 
physical  grounds.,  further  conditions  are  imposed. 

Thus,,  for  scatters  having  rotational  symmetry,, 

V(r^r*!)  can  only  depend  on  r  and  r 1 .,  and  on  the  relative  angle 
cl  between  r  and  r but  not  on  the  individual  orientation 
of  the  two  position  vectors.  A  physically  meaningful  kernel  can 
therefore  be  expanded  in  the  form 


(II.4) 
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where  the  factor  (2 f  +  1)/(4tt  rr 1 )  has  been  introduced  for 
convenience.  The  functions  Kg(r,r!)  determine  the  angular 
momentum  dependence  of  the  potential. 

Within  Quantum  Mechanics,,  nonlocal  (NL)  potentials 
have  been  used  mainly  in  the  study  of  nucleon-nucleus 
scatterings.  It  is  of  interest  to  see  how  they  formally  arise, 

(7) 

and  to  discover  some  of  their  basic  properties  (Feshbach,  1958) . 

We  shall  be  as  brief  as  possible. 

Consider  a  system  of  A  +  1  nucleons,  consisting  of  a 
nucleon  plus  a  target  nucleus  of  mass  number  A.  The  states 
of  the  target  nucleus  are  described  by  a  complete  set  of 
functions  r^),  corresponding  to  energies  .  For 

simplicity  we  neglect  the  spin  variables.  The  complete 
wave  function  f  can  then  be  expanded  as  follows : 


=  S 
i 


^  c 


i 


u 


(r  ) 
v  o ' 


(II. 5) 


where  r  refers  to  the  incoming  nucleon.  The  term  i  =  0 
describes  the  nucleus  in  its  ground  state  plus  a  nucleon 
in  the  incident  channel.  Let  denote  the  total 
Hamiltonian : 


HT  HA  ^ rl J 


r n )  +  T  +  V(r  , 
A'  o  v  o' 


where  is  the  Hamiltonian  for  the  A  particles  of  the  target 

nucleus,  Tq  the  kinetic  energy  operator  for  the  variable  r^, 
and  V  the  potential  of  the  nucleon  in  the  field  of  the  nucleons 
of  the  target  nucleus. 
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Substituting  expansion  (5)  into  the  Schrftdinger  equation 


Ht  V  =  E  f 


where  E  is  the  total  energy,  and  using  the  relation 


Hn^.  =  e.if/. 
Ari  iri 


as  well  as  the  orthonormal  properties  of  the  set  we 
obtain  a  set  of  coupled  equations : 


(T  +  V.  .  +  e.  -  E )u. 
v  o  11  1  '  1 


V.  .u  . , 

jFi  ij  j 


where 


V.  .  (r  )  =  <f.  |  V  |  ip  .> 
ij v  oy  i  r j 


To  derive  an  uncoupled  Schrddinger  equation  for  u  ,  we 
introduce  the  matrices 


u 


u  = 


u 


1 

2 


W 


(ii 


Vo  =  (Vol’  V'  •••  }  ’ 


(II 


/\ 


and  a  matrix  operator  K  defined  by 


/\ 

K  .  =T  *5..  +  V .  .  +  e  b.  . 
ij  oij  ij  iij 


(ii. 6) 


7) 


8) 


(II. 9) 
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Then  eq.  (6)  takes  the  form 


+  V 
oo 


/A  /\ 

V  u , 
o  J 


(TI.lQa) 


(H.lO.b) 


A+  A 

where  denotes  the  Hermitian  conjugate  of  V  .  Solving 

eq.  (lO.b)  formally  we  obtain 


where 


y\ 

u 


E^  -  K 


u 

o 


=  E  +  IT) 


(II. 11) 


(11.12) 


The  term  lr\  ensures  that  only  outgoing  waves  will  be 
present  in  the  exit  channels  u^  (i  >  l).  Eliminating  u 
from  eq.  (lO.a)  by  means  of  eq.  (ll),  we  obtain 


+  V  +  V 
oo  o 


E 


-  E)uq  =  0  (11.13) 


The  effective  potential.,  u  ,  is 


x> 


V  +  V 
oo  o 


1 


T+T 


E 


/A 

-K 


(II. 14) 


Several  interesting  general  properties  of  u  emerge  from  a 

study  of  eq.  (l4)  .  First  o  is  complex,  due  to  the  term  irj 

(+ ) 

contained  in  E'  The  complex  term  arises,  physically. 
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from  processes  in  which  the  nucleon  leaves  the  entrance 

channel  u  and  is  emitted  in  one  of  the  exit  channels.,  u.  . 

o  1 

Of  course,,  for  this  to  happen,,  one  must  have  E  >  e.  so 
that  reactions  could  occur. 


Secondly,,  r  is  nonlocal  in  character.  To  show  this. 


/N 


let  us  insert  the  eigenfunctions  of  K  into  the  space  between 


the  operators 


1 


T+T 


,  ,  ,  /N 

E  v  '  -K 


and  vt  in  eq.  (l4).  In  general. 


these  eigenfunctions  will  consist  of  a  discrete  part  and  a 
continuum,  {0^}  and  corresponding  to  energies 

{6^}  and  {£},  respectively;  a  is  the  degeneracy  index. 


We  then  obtain 


u  =  V  +2 
oo  n 


V  0  >><<£  V ^ 
on  no 


r 


E  -  £ 


+ 


n 


J 


da 

J 


V  ®(£,a)X$(£,a)v1 

at  -5 - m - °- 

E(+)  -  £ 


(11.15) 


Consider,  for  example,  the  second  term.  Writing 


the  matrix  element  between  I  r  >  and  |  u  )  ,  of  the  second 


term 

above. 

becomes 

<r 

v  0 

.  a 

V  0 
o  n 

><0  vt  |  u  > 
n  o  o 

n 

E 

-  e 

n 

—  y 

y 

f  V  .  ( r  )b.^n^(r  ) 
01 v  O '  1  v  o ' 

[h(n)T0')  5  VOU0(O  ’ 

Zj 

n 

Zj 

i ,  jVo  J 

C  0 

E  -£n 

;  ■  i  *  '  '  •:  ■  ■ 
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which  is  of  the  form 


r 


L(r  .  r  ' )  u  (r  ’)  dr  1  , 
y  o  o'  o  '  o  '  o 


i.e.,  nonlocal  except  in  extreme  circumstances 


We  note  that  the  numerators  in  the  expansions  of 
eq.  (15)  are  positive  definite.  We  may  decompose  o  into 


the  form 


0  =  0.  +  i 

1  2 


(11.16) 


where  and  are  given  by 

t 


V  $  ><0  V 

0  =  V  +  Z  - 8_°_  +  P 

1  00 


d  e 


r 


E  -  £ 


E-£ 


n 


da  VoO(£,a)><0(f  ,a)  Vq 


r 


and 


=  ~7r 


da  V  ,  a)><<E>(E,  a)  V 


t 


E  >  e 


1 


=  0 


E  <  e 


(II.18) 


1 


In  eq.  (17)  the  symbol  P  stands  for  Cauchy  principal  value. 


(11.17) 


We  note  that  is  a  monotonically  decreasing  function  of  E: 


80. 


W~  - 


<  0 


(11.19) 


and  that  is  negative  definite  --  as  would  be  expected, 
since  physically  this  means  that  the  incoming  nucleon  density 
may  be  decreased  but  not  increased. 
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Another  important  point  which  is  frequently  missed 
should  also  be  discussed.  Eqs .  (17,18)  show  that  u  and  o2 
are  Hermitian.  If  now  we  take  the  Hermitian  conjugates  of 
both  sides  of  eq.  (16),  we  find 

^  =  v  -  iu2  .  (11.20) 

Thus  0  is  Hermitian  if  and  only  if  =  0  identically. 

The  Hermitiari  condition 

-*■  1  ^  * 

<r  I  u  I  r'>  =  <r'  I  u  |  r  > 

is  therefore  unduly  restrictive,  and  in  fact  incompatible 
with  the  optical  model. 

II. 3  Symmetry  of  Kernel  function. 

( 10 } 

Within  certain  general  formulations v  of  the 
nucleon-nucleus  scattering  problem,  which  encompass  the 
highly  successful  optical  model,  one  finds  that  it  is  best 
to  deal  with  symmetric  potentials.  The  configuration  space 
is  divided  into  two  regions.  The  first  is  an  internal 
region  in  which  all  the  nucleons  interact  strongly  and  form 
nearly  stable  configurations.  The  second  (external)  region 
consists  of  channels  corresponding  to  the  various  ways  in 
which  the  compound  nucleus  can  disintegrate.  The  dividing 
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surface  is  roughly  a  nuclear  surface.  A  complete  set  of 
wave  functions  for  the  internal  region  are  obtained  by 
insisting  that  they  satisfy  the  many-body  Schrddinger  equation 
and  homogeneous  boundary  conditions  at  the  dividing  surface 
r  =  R.  The  wave  function  describing  the  actual  reaction 
process  may  then  be  described  in  the  internal  region  by  an 
expansion  in  this  complete  set,  and  in  the  external  region 
by  the  known  wave  functions  for  the  disintegrated  system  in 
each  channel.  Thus  in  the  internal  region  we  need  a  set  of 
functions  whose  radial  parts  are  solutions  of  the  eigenvalue 
equation  [see  eq.  (111,3)]. 


dr 


]  u(m)  (r) 


r 


For  simplicity  let  us  consider  the  S-wave  equation 


dr 


(r,  r '  )  u«  (r  1 )  dr 


(11.22) 


where  u 


According  to  the  formulations  of  the  problem  (for 
definiteness,  that  of  Kapur  and  Peierls^1^)  the  boundary 


conditions  are 


fm)(0)  =  0 


(11.23) 


[r  u(m)]  „  =  ikRif (R)  , 

dr  r=R  '  7 


> 


J 


’ 


fa  Iv  15  91 


'  J  :'nr  ..  nc  r  .  •  •' 
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and  the 
sense  : 


are  required  to  be  orthogonal  in  the  following 


(m)  (n) 
u v  'u v  ' 


dr  =  0 


m 


^  n 


(11.24) 


(n ) 

Multiplying  eq.  (22)  by  uv  '(r),  integrate  and  then  subtract 
the  corresponding  equation  obtained  by  interchanging  m  and  n, 
we  obtain  with  the  help  of  boundary  conditions  (23) : 


(E(m)-E(n))  /fn)UW  dr  = 

//  dr  dr'  Kq  ( r,  r '  )  [fn4  r,)i4m4r'  )  -iln4r '  )u  ^m4r)  3 

/  \ 

A  barely  sufficient  condition  for  the  orthogonality  of  u^  ’ 
and  u^n^  when  n  /  m  is 

K0(r,r' )  =  Ko(r',r)  , 

or,  in  general, 

K^(r,r'  )  =  K^rqr).  (11.25) 

From  eq.  (4)  we  deduce  the  equivalent  relation  to  eq.  (25) : 

V(r,rT)  =  V(rT,r)  (II. 26) 

The  above  argument  leading  to  eq.  (26)  is  essentially  that 
of  Brown  and  De  Domini cis  ^  . 

The  condition  of  symmetry  on  V(r,rT)  is  met  if 
the  matrix  elements  <r  |  u  |  rT>  and  <r  |  I  r»>  [eq.  ( l6 )  ] 


are  real . 


r'ri/j  c  o 

m)i  lc  ■  J  '  g<  -  to  .-irW  'Xol  o.  •  ;  cieiJ>xr'*J8  9;i£d  * 

jj  on..:, 

*  =  ( ,'iV*)03t 

tl£rit  n  nJt  •('to 

;  .  s  ,)  9  1  (  JisXSl  ft  9l  >ViUp9  9  I'D  >  OV  (•')  .  p'3  !®v;ri 

(ds.n)  (*,’*)v  = 

3  9ti  si  ('  <  n } V  no  Y‘1  arrriYS  lc-  it©J  -tbrceo  nT 
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We  note  that  for  a  complex  local  potential*  eq.  (26) 
is  automatically  satisfied.  Henceforth  our  attention  will 
be  confined  to  symmetric  potentials. 

II. 4.  Flux  Density. 

Multiplying  the  Schrddinger  equation  (2)  by  , 

and  then  subtracting  from  the  result  the  corresponding  complex 
conjugate*  we  obtain 

S'. 

V  •[  (W*  -  ****)]  =  J--  Im  /dr*  V(  r*  r T  )  f  ( r '  )^*(r)  (11.27) 

The  quantity  inside  the  square  bracket  is  the  usual  expression 
for  probability  current  density  (flux).  Noting  that  the 
right  hand  quantity  is  in  general  different  from  zero*  even 
for  real  V(r*"r')*  one  is  tempted  to  modify  the  mathematical 
expression  for  the  flux*  so  as  to  preserve  the  form 

V-T  =  0  (11.28) 

for  real  potentials*  where  J  is  the  current  density  vector. 

We  recall  that*  for  the  case  of  a  charged  particle 
moving  in  a  magnetic  field*  the  expression  for  flux  has  been 
modified  to  include  terms  involving  the  vector  potential* 
though  for  different  reasons . 


(dS)  . p9  XL .  i-n^ooq  l£!>cl  xelqmoo  £  *iol  scrorc  si* 

niB^cfo  ow  i9^63ut,noo 


(*?)**( -~)  V('~,*E)V  '•■■{■  i.  ml  -  {(^Hf  -  * i 7>Sf>)  £-  ].  V 


-x  )  £  r?j.r  e:  L  ,  s  I :  £•  d  9fi£jj]  en  t .!'  '  an  I.  v  ;  .;  c  d rfT 

,  rme‘t  t lb  .£•  erj^;  al  V  JE^-r >  L  £i 

03.60  SfCtf  T;Or  tJ.6rl^  Xj£0  3^t  sW 
XU  o']  ,3*  •  ;»jXf  Mi  ,  bl£  'r‘l  \Z  [£•  ft  "••fl-tVO 

t  i.s  -nsjc'!  ic  5  *,  O;  \1  jv  y\  ovr  rrrisl  <  Coni  oj  bs  ;  .  bom 
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On  the  other  hand,  when  the  potential  is  complex, 
actual  absorption  occurs,  as  in  the  case  of  local  complex 
potentials.  We  must  therefore  allow  for  the  appearance  of  sinks. 

The  common  practice  is  to  leave  the  expression  for 
the  flux  unchanged,  so  that  eq.  (27)  is  written 


(11.29) 


For  real  potentials,  appearance  is  saved  by  introducing  the 


interpretation,  flux  sometimes  leaves  the  elastic  channel, 
enters  the  invisible  channels,  and  reenters  the  elastic 
channel.  The  fact  that  over  all  flux  is  conserved  [as  can 
be  seen  from  eq.  (29)  since  V(r, r ’ )  is  symmetric]  is  because 
as  much  flux  departs  for  the  invisible  channels  as  returns 
from  them. 

It  appears  that  an  alternative,  which  satisfies 
eq.  (28)  when  V(r,r ’ )  is  real,  and  gives  actual  absorption 
when  V(r,r T )  is  complex,  is  possible.  One  has  only  to  modify 
"J  in  eq.  (29)  so  as  to  absorb  the  effect  of  the  real  part 
of  V(r,r ' ) .  We  thus  obtain  in  place  of  eq.  (29)  the 
following  equation  of  continuity: 


V  * 


1  =  J-  Re  /V2( ?,?')  V'(r')  ip*(r)  dr' 


(11.30) 


where 


n rc. \  w  :  ,  'v 

.. 

•e  r.-'  o . •  :j 

■ 

■TV1  )  V  ion  (  .  "  •  :*n  .  i'  9.: 

■  ..  l:  ;  v  1  -  9  :r  .  ■-  :  i  •, 1  "  ;  .  > :  j  n  r;« 

A-'  •  8  l  V.  trfw  dv.f^fi04l  ,cl£  f.  S  lB0r  -JS  *1 

r  -  q  J  ‘  r-  ,  >  v  r  .  'iog-:  J  c  3  .  >i  ru  \ 

(01. 1  )  (*)*",-  ( 4).*  (  ^*t)  ,V\  9S  :l-t  .  '  7 


i6 


7  =  -  ^  -  JU-  2  7/  dx^m/dr’  V]L  (r,  rT  )^(rT  )^*(r) , 

IL  1 

n  =  dimension  of  the  problem  =1,  2,  or  3, 
r  =  (x1,  xg,  x^) 

and  and  are  the  real  and  Imaginary  parts  of  V, 
respectively. 

Since  flux  is  not  observable  in  principle,  there 
is  some  arbitrariness  in  its  mathematical  definition.  In 
the  interest  of  simplicity  of  interpretation,  the  second 
alternative  may,  therefore,  be  tentatively  accepted. 

II. 5*  Wigner’s  Inequality. 

In  non-relativisti c  physics,  the  causality  condition 
states  loosely  that  a  cause  must  always  precede  its  effect. 

Thus,  in  scattering  problems,  the  scattered  waves  cannot 
appear  before  the  incident  wave  reaches  the  scatterer.  As 
is  well  known,  the  condition  leads  to  a  mathematical  inequality 
which  sets  for  each  scattering  problem  a  lower  limit  for 
the  rate  of  change  of  phase  shift  with  respect  to  incident 
energy  (Wigner’s  Inequality v  J  )•  It  might  appear  surprising 
that  the  causality  condition  as  stated  above,  being  in  fact 
a  tautological  statement,  could  have  a  consequence  of  such 


I ,  • 


. 

^  ,  0  tor?  xq  n  r  s  .cfBvneado  .“or  xu n  an.  F 
r  •  '  '•  ‘'3rf?£rn  ad.r  rri  e  >n  i..  i  .  -3J; 

b'  0  )P‘  ?JfW  'll<  *o  YjZsli  r;ils  J.o  se*  9.:  ii  ecit 


' 

*  05r:  3  2;:  1  -  ^eo  nq  3ybwJ\£  iaum  earso  b  xBnY  ^  a  >oo..  asdBda 
orrn.80  3  asw  Iwxetfjjoa  sdd  tem9ldo*q  ?.. xlnod. ..--a  ni  .  wri. 
3A  •  s‘  3':  'SZ'S  9ri:  a  mojBS'i  e-'Bv;  dTteJb  :onl  arid  enoAs  j  ojBeqqB 
/p'  -  Ib  tr  mtr  "i  Bffl  b  od  a6B9j  no  id:  -000  arid  ,,  ,;on:i  CXsv  g; 

d0Br  n  3fti9d  vSVoc <8  bacteda  as  noMibnoo  VdllBPJJBO  9rid  dBrid 
03  9  ean.^o  b  s  /b.  :  b.uoo  <  namedBda  Iroiso  odu.  d  b 
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physical  interest.  However,,  the  causality  condition  is  not 
so  much  a  condition  imposed  from  without  as  a  reminder  of 
the  concepts  of  cause  and  effect  already  incorporated  in 
the  formulation  of  the  problem.  As  such  it  has  the  status 
of  a  mathematical  identity. 

The  WignerTs  inequality  has  been  derived  for  real 
energy-independent  local  potentials.  Let  the  range  of  the 
potential  be  b.  Then 


sv 

dk  dr 

1 

U  *  •« 

&  dkdr 

>  0 

Wigner ’ s 

— 

r=b 

Inequality 

(11.31) 


where  r  u^(k;r)  is  the  i-th  radial  wave.  For  1  =  0, 


u  (r) 
o v  ' 


r>b 


=  const  x  sin  (kr  +  5) 


where  6  is  the  S-wave  phase  shift.  Then  eq.  (31 )  gives 


d6  x  ^  ,  sin  2(5  +  kb)  x 
3k  -  '  b  +  2k  ^ 


(b  +  ) 


2k 


We  shall  now  show  that  the  same  relation  (31 ) 
holds  for  real  energy-independent  NL  potentials  which  vanish 
outside  certain  ranges.  The  following  proof,  first  noted 

( ill ) 

by  Razavyv  ’ ,  is  a  generalization  of  the  one  given  by 

( 15) 

Corinaldesi  and  Zienauv  ' . 


■ 


<  < 


.'■a  9/1  i.  ;ri1  -c  ...  re n  Xj.  n.;  31 


18 


We  have  the  equation  [see  eq.  (ill. 3)] 

u^(r)  +  (k2  -  —■ s|+I  ^  )  Ug(r)  =  /  K^(r,r’)u^(r')drh  r<b 

where  b  Is  the  range  of  the  potential.  Taking  the  complex 
conjugate  and  the  derivative  with  respect  to  k,  we  obtain., 
respectively., 

(r)  +  (k2  -  u^*( r)  =  K^(r,r’  )  u^r’Jdr' 

and 

u^(r)  +  (k2  -  d^(r)  +  2k  u^(r)  = 

b 

/  ^(r,r'  )  u^(rT  )  drf , 

where  the  dot  denotes  differentiation  with  respect  to  k. 

We  multiply  the  first  equation  by  u^(r)„  the  second  by  -u^(r)„ 
and  add  up  the  resulting  equations,,  obtaining 

H  y.  M 

u^(r)u^*(r)  -  u^(r)u^(r)  -  2ku^(r)  u^(r)  = 

=  /  K^(r,r')  u*(r')u^(r)dr’  -  /  (r, r T ( r ’ )u*( r) dr ’ 


or 

(u^u^  -  u^  u») 


-i  i 


b 

2k  u!  u„  +  /  drT  K  -  ( r,  r T  )  [u  Q  ( r  ’  )u  0  ( r )  -u  g  ( r  1  )u 

0  30 


■£  i 


i 


i 


(r)] 


•  ■  1 


■ 


'  c  i  •: :)  vrij  bn;,  kj 


bn,£ 


e(r  ;  v,cf  rollup©  cs'i.tl  9(1j  ^cXc*± J Xur*  ©w 


t  I  -  ;  ’1  l),  H  \ 
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If  now  we  integrate  both  sides  of  the  equation  with  respect 
to  r  from  0  to  b,  the  result  is 

Cui(r)u^*(r)  -  u*(r)u^(r) ]r=b  =  2k  /^  |  u^  |  2  dr  >  0 
since  K^(r,rT)  is  symmetric  and  u^(0)  =  0  for  all  k. 


.  . . 
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CHAPTER  III.  S-WAVE  SCATTERING  (WITH  NUMERICAL  RESULTS) 

The  analysis  of  proton-proton  scattering  experi¬ 
mental  data  is  one  of  the  primary  means  by  which  physicists  try 
to  understand  the  nature  of  nuclear  forces.  As  yet,,  no 
adequate  theory  has  been  found. 

In  this  chapter.,  two  attempts  at  fitting  proton- 
proton  "^"S  phase  shift  data,  for  the  energy  range  20  -  3^0  MeV, 
will  be  made.  The  potentials  used  are  (in  parts)  nonlocal. 

They  have  the  (rare)  merit  of  being  soluble,  at  any  rate  for 
the  S-wave  equation. 


III.l. 

Consider  a  system  of  two  particles  1  and  2  inter¬ 
acting  with  each  other.  The  total  Hamiltonian  operator  is 


+  /: 


] 

V 

R= 


where  V  characterizes  their  interaction;  the  other  symbols 
are  either  standard  or  self-explanatory.  For 


)] 


we  have  local  interactions. 


p| 


iO  a  -a'  (.'■  >;  I 


,0;  ..f'-.’i  n  ,il  ?j  :■ 
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The  motion  may  be  separated  into  two  parts  consist¬ 
ing  of  the  motion  of  the  center  of  mass  and  their  relative 
motion.  The  deduction  is  similar  to  that  of  the  case  of 
local  interactions.  The  Schrodinger  equation  for  the  rela¬ 
tive  motion  is 


(V2+k2)  f{v)  =  ^  f  V(r,r')  f(r')  dr’ 

ic 


(III  .1) 


where 


Tik  =  initial  relative  momentum 
and  m  is  the  reduced  mass. 


We  expand  the  wave  function  ^(r)and  the  kernel  function 
V(rqr')  in  terms  of  spherical  harmonics: 


f(r) 
V(r}r'  ) 


=  nZ  a 


Vr) 


m 


£7m  r  (9jtp) 


2£+1  Kj,(r,r'  ) 


Htt 


rr 


P^cos  a) 


=  „Z 


K£(r,r' ) 


Y,m(0.q>)  Y  m*( e  ' , cp  ’  ) 


£  ,m  rr T 


£ 


£ 


(III. 2) 
(II  u) 


Substituting  these  expansions  into  eq.  (l)  and  making  use 
of  the  property 

(0,,cp’)Yim(eT,cp,)dn'  = 


5  6 

££'  mm' 


n>  -  -  j  ■  •.  ■  ; 
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we  obtain 


Z  <  + 

£  ,m 


id  Q  d 

hits-  30  •  sin  e  3e 


+ 


1  d 


2 


sin20  dcp2 


Y/>^aO 


,  ^2  ^  m  aiu^  v  m  2m  r  TA  _,N 
+  k  Y„ - Y „  —  J  K  (r,  r  ’  ) 


£  r  i 


di' 


a^(r') 


dr 


=  0  . 


We  then  deduce,,  in  the  usual  way,  the  equation  for  the  1-th 
partial  wave: 


u£  +  (k2-^~J-  )  u^(r)  -  /  K^(r,r'  )  u^  ( r  T  )  dr T  =  0  (ill. 3) 


III. 2 

Consider  a  potential  of  range  b,  whose  partial 
kernels  K^(r,r' )  are  solutions  of  the  differential  equation 

2 

d  o  Kdr,r'  )  -  (p2+^(^.)  )  K,(r,rT)  =  -7\^26  ( r-r  ’  ) , 

dr2  1  r2  £  (III. 4) 

r,r'  <  b 


with  homogeneous  boundary  conditions 


h (O’ r ' )  =  K^(b,r')  =  0 


(III. 5) 


J  =  'ib  >1  \  |  -  (i)  «  ,  (■!),”« 


:  ov^'w  is sq 


< 
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This  example,  a  special  case  of  a  more  general  one 

( 14 ) 

considered  by  Razavy  ,  is  instructive.  Here  the  real 
constant  A  in  a  sense  (as  we  shall  see  presently)  measures 
the  strength  of  interaction,,  while  l/p  (>  0)  is  an  index  of 
"nonlocality "  in  the  sense  that,  as  l/p  0,  the  potential 
becomes  local.  For,  taking  the  limit  p  -*  00  ,  we  have  from 
eq.  (4)  that 

(r,  r  1  )  =  A5  ( r-r T  ) 
so  that  eq.  (3)  becomes 

Ug(r)  +  (k - 2“  )  up(r) - 2  Aup(r)  =  0  *  r  <  b, 

r  *fi 

which  is,  in  fact,  that  describing  a  potential  well  of  depth  -A. 

Since  potential  well  is  known  to  be  a  rough  approxi¬ 
mation  to  the  actual  interaction  between  protons,  it  is  worth¬ 
while  to  see  if  this  slightly  more  sophisticated  potential 
(plus  a  small  hard  core)  is  any  better.  The  hard  core  is 
added  to  give  short-range  repulsion. 

Let  us  first  solve  the  problem  without  the  hard  core. 
The  S-wave  equation  is,  from  eq.  (3) 

u"(r)  +  k^u(r)  =  f(r) 


r  <  b 


(III. 6) 


v \ j  .  i  ,d  fosis^lanoo 


; 


■tJ8rf,*t  ( &)  ,  p 9 


3  3r  o  o  ed  ( £ )  .  p  9  d  >3  fir  c  3 
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where 


and 


Oryi 

f(r)  =  -p  /  K  (r,rT)  u(rT)  drT  (ill. 7) 

Ui  ° 

u(r)  =  uQ(r)  • 


Multiplying  on  both  sides  of  eq.  (4)  (for  £=0)  by 

u(rT)  and  then  integrating  with  respect  to  rT ,  we  get 

f " ( r )  -  P2f(r)  =  -  %  Ap2u(r)  (ill. 8) 

-ir 

Eq.  (8)  can  be  used  to  eliminate  f(r)  from  eq.  (6).  The 
result  is  a  fourth-order  differential  equation  for  u(r) : 


4 

d  u  ,2  .  2v 
— -tf  -  (P  -k  )  2 

dr  dr 


,2 

d  u 


-  p2(k2-%  A)  u  =  0 
tfi 


r  <  b  (III. 9) 


The  boundary  conditions  (5)  are  equivalent  to 

f(0)  =  f ( b)  =  0  (III. 10) 


In  addition,,  we  require  the  usual  regularity  condition 


u (0 )  =  0  . 


(Ill  .11) 


Using  eq.  (6)  we  therefore  have  the  boundary  conditions 


u ( 0 )  =  u"  (0)  =  0 


2 


V. 


u"  (b-)  +  k  u(b-)  =  0 


(III .12) 


j 


’ 
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The  solution  of  eq.  (9)  is 


:u  =  A  sinh(pr)  +  B  cosh(qr)  +  C  sin(vr)  +  D  cos  (vr),  (ill. 13) 


where  the  constants  A,  ...  ,  D  are,  as  yet,  arbitrary,,  and 


li  =  — \  (p2-k2)  +  [(p^-k*)^  +  A)] 


2,2x2.  i.  „2  /,  2  2m 


1 

2 


■ft 


2 


i(HI.lA) 


V  ;  = 


1 


2  rt2 


(kc-p^)  +  [  (p2-k2)2  +  4p2(k2-  ^  A)] 

fi 


y/2 


We  shall  be  interested  in  negative  A,  so  that  the  overall 
effect  of  the  potential  is  attractive. 

Outside  the  range  of  interaction, 

u(r)  =  const  x  sin  (kr  +5)  r  >  b  (ill. 15) 


The  wave functions  in  the  two  regions  are  joined  smoothly: 


Ub-)  =  u(bf) 

u ' (b-)  u ' ( b+ ) 


(III .16) 


Three  of  the  constants  in  eq.  (13)  can  be  eliminated 
by  the  boundary  conditions  (12).  Then  eq.  (l6)  gives  the 
desired  relation  for  the  phase  shift  8  as  a  function  of  k: 


i ■  s  0  +  (>!i))ne<  9  •{  l  t)rii.j'F  A  y 


■  ...  -  ■  »rij  yiariw 


Vv 
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k  cot  [6+kb]  =  -g-1-  g  [(ii2+k2)v  cot ( vb)+  ( v2-k2) p  coth(ab)  ] 

li  +v 

(III. 17) 

If  a  hard  core  of  radius  c  <  b  is  added,  condition  ( 11 )  must 
be  replaced  by 

u ( r)  =  0 ,  r  <  c  (III .18) 

Let  us  also  modify  the  boundary  conditions  (5)  for  ,  to 
read 

K0(c,r'.)  =  Ko(b,r’)  =  0  (ill. 19) 

It  follows  that  the  boundary  conditions  for  u(r)  in  the 
region  c  <  r  <  b  are  obtained  from  (12)  by  replacing  u(0) 
and  u"(0)  by  u(c)  and  u"(c),  respectively.  We  obtain  an 
equation  slightly  different  from  eq.  (17),  viz. 


k  cot  (kb+b)  = 


1 


2  2 
|1  +v 


(p2+k2)v  cot  [  v  (b-c)  ]+  (v2-k2)a  coth  [pi(b-c)  ] 


(III. 20) 

With  eq.  (20),  we  try  to  fit  the  S  phase  shift  to 

Breit's  values  (YLAM)^1^^,  by  adjusting  the  various  parameters. 

We  find  that  the  following  set  of  values,  though  perhaps  not 

the  best  choice,  gives  a  reasonable  fit.  These  values  are 

(17) 

close  to  those  of  Giltinan  and  Thalerv  ' . 


7\  =  -27  MeV,  6  =  2.174  F  1,  b  =  2.6  F,  c  =  0.157  F. 

I* 


The  result  is  given  in  table  I,  under  the  potential  u 
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Elab(«eV) 

20 

100 

180 

260 

340 

6  (u^ )  (rad . ) 

0.884 

0.300 

0.038 

-0.090 

-0.157 

5  {vxl) (rad. ) 

0.857 

0.341 

0.094 

-0.071 

-0.195 

5 ( YLAM) (rad. ) 

0.856 

0.380 

0.136 

-0.033 

-0.195 

Table  I . 


III. 3. 

Consider  another  potential*  consisting  of  a 
nonlocal  one  of  range  b  given  by  eq.  (4),  and  the  Eckart 

potential  P  _ 

2  V  7  e  y 
o  ' 


"  '(1+V„fr)2 


for  r  >  b. 

We  shall  confine  the  parameters  (3  and  A  to  values 
satisfying  the  inequalities 


4e2(%  a 

-fi 


k2)  >  (p2 

P2  >  k2 


k2) 


(III .21) 


for  the  energy  range  20  -  340  MeV.  This  means*  in  particular* 


that  we  want  the  nonlocal  potential  to  be  repulsive. 


For  r  <  b*  the  differential  equation  is  identical 

with  eq.  (9)*  Its  solutions  are  linear  combination  of  exponen- 
or 

tials  e  where  a  are  roots  of  the  algebraic  equations 


.  9ii  ::T 
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2a2  =  (p2-k2)  ±  i  /4p2(2 ™  A  -  k2)  -  (p*-k*) 

-fi 


,2  ,  2xP 


The  aTs  are  given  by 


cx  —  p  +  i^,  -p  4-  i<^ 


(III .22a) 


where 


P  =  \  J  2p  J  -  k2  +  (p2 


0)  =  2p  -  k2  -  (p2 


k2) 


k2) 


2 


i  (III. 22b) 


Next  we  employ  the  boundary  conditions  (12)  to  eliminate 
three  of  the  four  constants  appearing  in  u(r) .  The  result  is 


u(r) 


_  _ sinh  (pr)  cos  (^>r) _ 

(p2-CD  +k  )  sinh (pb)  cos  (<^b)  -  2pu>  cosh(pb)  sin(<^>b) 


_ cosh  (pr)  sin  (<^>r) _ 

(p2-o2+k2)  cosh(pb)  sin(^b)  +  2pco  sinh(pb)  cos(^b) 


where  we  have  set  the  remaining  constant  equal  to  unity,  as 
we  are  mainly  interested  in  the  logarithmic  derivative.  The 
value  of  this  latter  quantity  at  r  =  b  is  then  found  to  be 


F 


u  T  ( b— 

u  ( b- 

P  (p^+^+k2)  fcoth(pb)  -  tanh  (pb)^+^>(p^+^>2-k2 


)  ('cot (^>b)+tan(^b 


]) 


2pco 


cot  (<^b) 
coth ( pb) 


+ 


coth 

cot 


(111.23) 
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For  r  >  b,  the  differential  equation  is 


2V72e)"7r 

u"(r)  +  k  u(r)  =  - 


,  -yr\2 

(1+  Vn  e*  y  ) 

This  has  solutions  of  the  forrr/ 


u(r) 


o 

(20) 


(III .24) 


u r>b  =  C  +  D  h(r) 


(111.25) 


where  C  and  D  are  arbitrary  constants , 


g  ( r )  =  jl  +  — 


1+4 


< 


k 
7 


h  ( r)  =  J  T  + 


Y(r) 


1+4 


k 

7 


sin  kr  + 


2  k/y 


1+4 


k 

7 


Y(r)  cos  kr 


g  +7 


1+4 


k 
7 


(III. 26) 


cos  kr  -  - - Y(r)  sin  kr 


and 


As  r 
form 


00 


Hence 


2  V  e"^r 
Y  ( r )  =  -  VQ 


1  +  V  e 
o 


r7r 


(111.27) 


,  Y(r)  ->  0  and  so  u(r)  has  the  required  asymptotic 


u(r)  ~  C  sin  kr  +  D  cos  kr. 


=  tan  & 


(III. 28) 
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The  logarithmic  derivative  is 


u *  ( b+ ) 
u(b+) 


C  g  T  (b)  +  D  h  1  ( b)  _  g ' ( b)  +  tan  5  •  h  T ( b ) 
C  g(b)  +  D  h  ( b)  ~  g  (b)  +  tan  b  •  h  (¥) 


The  logarithmic  derivative  must  be  continuous  at  r  =  b.  In 
terms  of  F  =  ^ as  defined  in  eq.  (23) ,  we  therefore  obtain 


or  else. 


tan 


r  _  S'  (b)  -  F  g(b) 
Fh ( b )  -  h'  (b) 


> 


tan  (kb  +  6) 


tan  kb  +  tan  5 
1  -  tan  kb* tan  5 


j 

[F  sin(kb)h(b)  -  cos (kb) *g(b) }  +  (g,(b)cos  kb-h'(b)  sin  kb] 

[F  cos(kb)h(b)  +  sin (kb) *g (b) }  -  [gf(b)sin  kb+hT(b)  cos  kb] 


Making  use  of  eq.  (26)  defining  h(r)  and  g(r),  we 
finally  obtain 


K  -  K  F 

tan  (kb  +  6)  =  — - - 

K2F  -  K4 


(111.29) 


where 


1  +  - - — 2“  Y(b) 

1  +  4 

7 

2  ^7  p  -  Y ( b) 
k 

1  +  *  ^ 


i 
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(hi. 30) 
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k 


1 

1  +  ,2 

k 

1  +  4  ^2  Y(b) 


Y '  (b) 

1  +  4 


-  2  |  Y(b) 


Y 


2  k/7 

7 


Y  (b)  ' 


(HI. 30) 


and  Y(r)  is  as  defined  in  eq.  (27). 

With  eq.  (29)  we  try  to  fit  the  lS  phase  shifts 
to  Breit’s  values  (YLAM),  by  adjusting  the  various  parameters , 
subject  to  the  restriction  (21).  We  find  that  the  following 
set  of  values  gives  a  reasonable  fit: 

A  =  788  MeV,  p  =  4.3  F'1  ,  b  =  0.7  F 
V0  =  73  Mev,  7  =  2  P'1  . 

The  result  is  given  in  table  I,  under  the 
II* 


potential  u 
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CHAPTER  IV.  EQUIVALENT  POTENTIALS 


IV. 1  Momentum-dependent  Potential. 

Nonlocal  potentials  are  mathematically  equivalent 
to  momentum-dependent  potentials'1^.  That  is  to  say,  for  any 
NL  potential  operator  V  characterized  by  its  matrix  elements 
V(r,rT)  =  <r  |  V  |  ?T>,  a  function  U(r,p)  of  the  position 
and  momentum  operators,  r  and  p,  can  be  found  such  that, 
formally,  the  following  holds : 


U(r,p)  ^(r)  =  /v(r,rT)  ip  (r T  )  drT 


(IV. 1) 


For  we  can  expand  the  wave  function  ip(  r  ’  )  appearing  in  the 
integral  about  the  point  r '  =  r,  thus. 


ip  ( r1)  =  ip(  r  +  r’  -  r) 


=  Z  ^-r  [  (r  ’  -  r)  •  V]  ip( r) 


n 


n  J 


=  2  C(r}  -  r)*i£  ]n  ^(r) 


=e 


-fi 


( r 1 -r ) * p 


ip(r) 


Therefore 


U(rr,p)  =  fdr*  V(rr,r’)  e 


i_ 

- fi 


( r ’ -r) • p 


(IV. 2) 


As  an  example,  consider  a  potential  V^,  having  Majorana 
exchange  character,  defined  by 


A  I  T&TqAHD 
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<?  |  VM  |  =  VM(r)  f(-r) 

Since  <r  [  VM  |  f  >  =  /dr '  <r  |  VM  |  r’>  <rT  |  , 

it  follows  that  <r  |  |  r’>  =  V^r)  8(r  +  rT ) 

Using  eq.  (2),  we  find  the  equivalent  momentum- 
endent  potential  to  be 

->  ~h  2r*p” 

Um(Up)  =  VM(r)  e 
The  example  is  due  to  Wheeler^1! 

IV. 2.  Local  Static  Potentials  in  Born  Approximation. 

Equivalent  field-theoretical  nuclear  potentials  are 

defined  as  those  which  give  the  same  K-matrix  elements  in 

( 19) 

the  first  Born  approximation v  ' .  In  the  non-relativisti c 
reduction  they  may  be  taken  over  and  used  in  the  Schr<3dinger 
equation . 

Given  a  NL  potential  whose  S-wave  kernel  K(r,r’) 

(eq.  III. 3)  with  &  -  0)  is  energy-independent,  we  wish  here 
to  find  a  local  static  potential  which  gives  the  same  S-phase 
shift  in  Born  approximation. 

The  phase  shifts,  in  Born  approximation,  are  given  by 

00 

tan  5  =  -  — /  dr  sin  kr  “o(r) 

"n  k  0 


(iv. 3) 
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and 


00 


tan  8 


NL 


2m 

p  j  VAJ- 1  ' 
-fi  k  0  1 


00 


/  dr.,  sin  kr^f  K(rn  ,,  rQ)  sin  krQdr 


0 


12 


2  2 


(iv. 4) 


Eq.  (4)  is  easily  derived  by  converting  eq.  (ill. 3)  for 
l  =  0  into  the  following  integral  equation: 


2m 


00 


00 


u(r)  =  sin  kr  -  — p- 

41  0 


f  G(r,r')  f  K(r 5 , r" )u(r" )dr !dr" 

0 


where 


1 


G(r^r' )  =  £  sin  kr^  cos  kr  ^ 


Here  r^  denotes  the  lesser  of  r  and  r  ! .  The  first 


Born  approximation  is 


u(r) B  =  sin  kr  - 


2m 

-R2 


00 


00 


/  G(r,r!)  /  K(rT ,r' )sin  krTdrTdr"  . 


0 


0 


Its  asymptotic  expression 

00  00 

u(r)^~sin  kr  -  cos  kr  *  *  i-  /  dr'sin  kr/  K(rT,r")sin  kr"dr" 

B  hi  k  0  0 


implies  eq.  (4) . 


We  ask  the  question:  given  that 

00  00 

f°°  dr,  sin^kr.  v>  ( r,  )  =  /  dr  sin  kr,  /  K(r,  ,rQ)  sin  krQdrp  (IV.  5) 
01  11  O1  1  0  1 

for  all  k  (and  hence  for  large  k) ,  what  is  the  functional  f(K) 

where  f(K)  =  x>  ? 
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First  we  note  that  condition  (5)  is  equivalent  to 


£  00  2  $ 

dr]L  sin  kr^r^)  =  ^  f  dr^sin  k^  /  K(r1,r2)sin  krQdr, 


00 


00 


2“x2 


(IV. 6) 

for  all  k  since  both  sides  of  equation  (5)  vanish  for  k  =  0. 

The  left  hand  side  of  eq.  (6)  gives 


00 


Ik  /  dr  sin2  kr  u(r) 


00 


f  dr*2  sin  kr'cos  kr*r  u(r) 
0 


0° 

=  :J  q  dr'sin  2  kr  r  'o(r) 


i-  /  dr* sin  kr  .r  *o (■£■) 
0 


00 


=  sine  transform  of  r  u(^) 

The  left  hand  side  of  eq.  (6)  gives,  after  some  deductions. 


00 


00 


BE  f0  drlSln  krl  /  K(rl'r2>  sln  kr2  dr2 
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00 


=  sine  transform  of  ^-r  /  drT  K(r-r’,rT)-r  drT  K(ixfrf,rT) 


0 


Hence,  taking  the  inverse  sine  transforms,  we  obtain 
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almost  everywhere,  or 

2r  oo 

u(r)  =2  /  dr’  K  (2r-rT,  r’)  -4  /  dr1  K(2rH-rT,rT) 
0  0 


As  an  example,  let 

K(r',r") 


-(r'+r") 

e 


Using  eq.  (7),  we  get 

u(r)  =  2 e  "2r(l  -  2r) 

Equation  (5)  is  easily  verified. 


(IV. 7) 


iv. 3. 

In  a  situation  where  the  fact  that  p”  and  r  do  not 
commute  may  be  neglected,  thereby  reducing  the  description 
to  a  classical  one,  it  is  natural  to  ask  for  the  (in  some 
sense)  classical  equivalent  of  the  potential  operator  V. 

We  may  take  TJ (r, p")  as  given  in  eq.  (2)  to  be  such  a  one, 
treating  p  in  the  integrand  as  a  C-number  variable;  or  else, 
following  Dirac we  may  again  choose  U(r,p)  but  with  r 
defined  as  the  mean  of  the  position  variables  appearing  in 
the  integral,  thus: 

Ucl  =  /  V(rT,rM)e“1^r  "r  ^  P/^d(  r 1 -r"  )  (IV.8) 


where 
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r  =  75-  (?'  +  r " ) 

Let  us  consider  a  potential  whose  partial  kernel 
functions  K^,(r,rl)  (eq.  II. 4)  are  solutions  of  the  differ¬ 
ential  equation 

^~2  tK^(r,r')  V(r)  ]  -  (p2/fi2+  Vi+Dj  L  (r,r')o 


-  -A  (32/f)2o2  (r1)  B(r-r') 


(iv. 9) 


with  boundary  conditions 

Kpo,r')  =  K^(«,r')  =  0  (IV. 10) 

This  potential  has  the  property  that,  as  f3  ->  00  ,  it  reduces 

to  a  local  potential  Au(r).  The  solution  of  equ.  (9)  is 

^  ^  ,  (14) 

found  to  be 

Kjg  ( r J  r r rr’  u2(r)  u2(r»)  j  ^  ( i^r<  )Y\  r>  ) 

hi 

(IV. 11) 

( 1 ) 

where  and  h^y  ’  are  spherical  Bessel  and  Hankel  functions 
of  the  first  kind,  respectively. 


Then 
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V(r,r')  =  -  rr'-K3  AP3  u2(r)u2(r')  2  (2^+1 )  j^l^ry) 


hf'  (if  r>)P^(cos  «) 
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(IV. 10) 


Let  us  set  v(r) 


e  ->r/a 


Prom  eq .  (8) , 


UC1(?,P)  =  /  V(r’,rTI)  e  "^r  r  )  P/^d  ( r  ’ -r " ) 
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Since  the  factor 
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fi  ensures  that  the  major 


r  -r 


contribution  to  the  integral  comes  from  the  region  in  which 
rtC'  r"  ,  the  effect  of  replacing  r’  +  r"  by  |  r’+r"  |  should  be 
small.  This  is  done  to  obtain  a  closed  form,  thus  making  the 
result  physically  more  transparent.  Then 
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uo1(?,p) 


We  observe 


A  -r/a  r 

w  ^  e  j 


e 


-  £  p 
-n  p 


e-ip-p/K  d?  ^ 


p  =  r  r  _rM 


-X  -  £  p  tt  -ip  £  cos  @ 

e  '  /  dp*p  e  2t r/  e  sin@  d0 

dn  0 


=  A 


'fip 


-r/a  n 
J  dp  e 

0 


P 


sinC^p) 


2  e  ~r/a 
A  p  -| - g 

P  +  P 


that,  as  p  -#■  oo 


Jcl  ->  A  e  r//a  =  A-o(r), 


as  should  be  expected. 
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CHAPTER  V.  KLEIN-GORDON  EQUATION 


V.l 

¥e  shall  take  as  our  starting  point  the  Klein- 
Gordon  equation  in  a  local  field  V(r) : 

(V2-m2)^(r)  +  (E  -  V(r)  )2  ^(r)  =  0  (V.l) 

where  E  is  the  relativistic  energy  of  the  one-particle  system. 
It  will  be  taken  for  granted  that  eq.  (l)  represents  a  kine- 
matical  improvement  over  the  usual  Schrftdinger  equation,  as 
a  description  of  high-energy  scattering  processes.  A  dis¬ 
cussion  of  the  many  difficulties  inherent  in  eq.  (l)  lies 

(21 ) 

outside  the  scope  of  the  present  workv  ' .  Our  object  is 
to  extend  eq.  (l)  to  include  nonlocal  interactions,  and  from 
there  to  derive  some  relations  of  interest. 

The  transition  to  the  case  of  nonlocal  inter¬ 
actions  is  equivalent  to  adopting  a  momentum-dependent  poten¬ 
tial  U(r,  p")  in  place  of  V(r) ,  since  formally 

U(r, -iV)^(r)  =  /  V(r,rI)  ^(r')  drr  ,  (V.2) 

or,  more  generally, 

U(r, -iV) f (r)  =  f  V(r,r')  f(r')dr'} 


(v.3) 


G  - 
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where  the  function  v(r,r'  )  is  determined  by  U(r,p),  and  f(r) 
is  a  function  of  r,  (See  Chapter  IV) . 

Then  eq.  (l)  becomes 

(V2-m2)  i/y(  r)  +  (E  -  U(r, -IV))2  ip{v)  =  0 
or  else,  making  use  of  eq.  (3), 

V2^(r)  +  (E2-m2)  f{v)  =  f  dr’  L(r, r’)  f(r')  (V.4) 

where 

L(r,r’)  =  2E  V(r,r'  )  -  f  dr"  V(r,r")  V(r",r’)  (V.5) 

We  note  that  L(r,r’)  is  symmetric,  since  V(r,r’) 
is  assumed  to  be  so.  (Chapter  II). 

Eq.  (4)  may  be  rewritten  in  the  form  of  an  integral 
equation.  Let  It  denote  the  incident  wave  vector, 

- - 

/  2  2 

k  =  J  E  -  m  , 

and  let 

p(r)  =  /  dr'  L(r,r’)  r '  )  (V.6) 

Then  eq.  (4)  becomes 

(V2  +  k2)  f  (r)  =  F ( r )  .  (V.7) 


Since  the  Green’s  function. 


.  3  c'’")-  ••  •  i  ■  ■  ;  "  * ) V  no  j  .\s,  f..«3  erroriw 
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G( v, r' )  =  - 


1 

Ijvr 


lk  r-rT 


r-r T 


(v.8) 


satisfies  the  equation 


(V2  +  k2)  G(r,r')  =  &(r-r!)  , 


(v.9) 


it  follows  immediately  that 


f{r)  =  ^  (r)  +  /  G(r,rT)  Fir1)  drT 

=  fQ(r)  -  ^  /  drT  - — -  f  dr”  L (r  ’  ,r")f(r") 

|  r-r!  | 

(V.10) 

is  a  solution  of  eq.  (4),  where  ^  (r)  is  any  solution  of  the 
associated  homogeneous  equation : 

*0(r)  =  C  eli?-?+  C2  e’1*'7  . 


For  scattering  problems,  we  take 
so  as  to  obtain  the  correct  asymptotic  form. 

f(r)  =  elIf'r  -  2 ^ - 1  /  d?" 

'  I  r-r’  | 


=  1, 
Thus 

L(r» 


c2  =  0, 

r" )  f(r"), 
(V.ll) 


The  integral  equation  (ll)  may  be  solved  by  iteration. 
The  first  Born  approximation  is 


. 


.  J! ' 
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f  (?)  =  e 


ik  •  r  1 


Wir 


I  dr 


ik  |  r-r  ’  | 


c  j  — T  /  — ► .  — >n  \  i  k  •  r 
J  dr  L( r ’  ,  r  )e 


r-r ’ 


=  elk’r  -  | —  f  dr ’dr" 


ik  I  r-r 


V(  r ’ ,r")  e 


ilc.r" 


r-r  ’ 


+  w 


f  drT  dr"  dr 


1 


ik  r-r 


r-r  ’ 


V(r',?1)V(r1,r")elk' 


(V.12) 

The  Born  approximation  has  been  tested  for  complex 
potential  well  (eq.  (hi  .4)  with  A  complex)  at  high  incident 
energies  (~  1  BeV),  by  calculating  both  the  exact  and  approxi¬ 
mate  S~wave  scattering  cross-sections.  It  was  found  to  be 
remarkably  good,  the  error  being  within 


v.3. 


The  Klein-Gordon  equation  (4)  can  easily  be  decom¬ 
posed  into  partial -wave  equations.,  in  a  way  completely  analO' 


gous  to  that  of  section  III.l.  Thus  writing 


f(r)  =  2  a  »— 

i.,m  r 


Ug(r) 


Y/  (e,v) 


„  2i+l  K4r’r'>  „  , 

V(r>r')=  2  -TrT  -  Pji  (cos  a) 


(V.13) 


(V.l4) 
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whence 

L(r,r’  ) 

=  ^irL/r4)pi  <cos  “)  (v-15> 

where 

L^(r,rT)  =  2E  K^(r,rf)  -  /  dr”  K^(r  ,  r"  )K^  ( r 1 ,  r”  )  (V.l6) 

we  obtain  for  the  i-th  partial  wave  the  equation 

[u^  +  k2  -  =  f  (r,r' )  u^(rr)  dr? .  ( V. 17 ) 


'  ■'  r  •  ■  C'.’fW  ■ 
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CHAPTER  VI.  SECOND  QUANTIZATION 

The  general  procedure  for  the  quantization  of  the 
equations  of  motion  of  a  classical  system  may  be  applied  to 
a  wave  field.  In  this  way,  for  example ,  one  obtains  the 
theory  of  photons  from  the  classical  electromagnetic  field. 

The  quantization  of  a  ^-field  such  as  that  described 
by  the  Schr<3dinger  equation  is  commonly  known  as  second  quan¬ 
tization.  The  transition  from  classical  particle  mechanics 
to  wave  mechanics  constitutes  the  first  quantization.  It 

/  -1  o  \ 

has  been  shown ^  '  that,  for  the  case  of  local  interactions, 

the  second  quantization  of  a  Schr<5dinger  ^-field  converts  a 
one-particle  theory  to  a  many-parti cle  theory.  We  shall 
generalize  the  argument  to  include  nonlocal  interactions  as 
well.  For  simplicity,  we  treat  the  one -dimensional  case. 

The  presentation  closely  follows  that  of  Schiff^1^. 

VI. 1  First,  we  postulate  the  existence  of  a  Lagrangian 
density  of  the  form 

oC  =  itf  f*f  ~  S  lx  !x“  "  /  K(x,y)^(y)dy  (VI.l) 

-00 

where  dot  denotes  differentiation  with  respect  to  time,  and 
K(x,y)  is  the  kernel  of  the  nonlocal  potential. 


x.  X*  >  '■  ).  0  -  -1  3j  ,oi,,  u  '  ' 
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The  Action  A  of  the  system  between  times  t^  and 
t£  is  defined  in  term,  of  the  Lagrangian  density  as 


t 


2  oo 


A  =  /  f  X 


dx  dt 


(VI. 2) 


t. 


We  vary  the  Action  A  with  respect  to  ip  and  ip  indepen¬ 
dently.,  subject  to  the  usual  restrictions 


BA(t1)  =  6A(t2)  =  0  ,  (VI. 3) 

and  employ  the  Principle  of  Least  Action  to  derive  the 
equations  of  motion. 

In  varying  with  respect  to  ip*,  the  functional 
dependence  of  o C  is 

£  =  £{  **,  ^  )  . 
and  so  the  Euler’s  equation  would  result: 


eq. 


d  _  c)  hj2  _  n 

W  '  3(^|)  "  3 

Substituting  the  expression  for  f (eq. 
(4)^  we  obtain 


(VI. 4) 


l)  into 


0° 

inify  -  f  K(x,y)  ip(y)  dy 
-00 


d  /  _  <h^  hip  \ 
£x  '  2m  c)x 
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or 


ttip  =  — |  -f-  f  K(x, y )  ^(y)  dy 

dx 


(VI. 5) 


which  is  the  time -dependent  Schrodinger  equation. 

In  varying  with  respect  to  ip,  we  first  note  that 
the  Lagrangian  density 

2 

X '  =  W  -  ~  T'W  /  K(x>y)  ^*(y)  dy 

which  can  be  obtained  from  oC  by  interchanging  ip  and  ip* 
in  the  last  term,  gives  the  same  action  as  does  JZ  ,  since 
K(x,y)  must  be  symmetric.  Hence  again  we  can  use  the 
Euler’s  equation  directly.  We  obtain  an  equation  conjugate 
to  eq .  ( 5) . 

2  2 

-i'Tty'*  =  -  —  +  /  K(x,y)  ^*(y)  dy  .  (VI.6) 

2m  dx 

The  same  equations  (5*6)  will  be  obtained  if,  instead  of 
ip  and  if/*,  one  treats  ip^  and  ip2  as  the  independent  fields, 
where  ip^  and  ip 2  are  defined  by 


ip  =  f1  +  lf2 
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VI. 2  The  Hamiltonian  density  of  the  system  is  defined  to  be 


where 


Hence 

U 


a,t) 

*  (  * 

=  7 Tip  +7 Tx  J  Ip* 

-/ 

(VI. 7) 

7 T 

^  ~  =  ifty* 

(vi. 8) 

S^ 

(  *) 
TTk  ' 

-  3f  -  0  . 

(VI. 9) 

9 

cty* 


2 

i +  0  -  +  755;  +  ^*(x)/K(x,y)^(y)dy 

lx  +  ^*(x)/K(x>y)  V'(y)  dy  (VI. 10) 


(  *) 

Note  that  Try  '  and  cannot  be  regarded  as  a  pair  of  canon¬ 
ically  conjugate  variables  since  they  cannot  be  made  to 
satisfy  the  usual  Poisson  equation.  This  gives  no  problem 
since  they  can  be  eliminated  from  the  Hamiltonian  density 
altogether.  Thus,  using  eqs,  (8,9)*  one  obtains  from  eq.  (10) 


2m 


St r 

c>x 


S 


i 

"tt  w 


/  K(x,y)  ^ (y )  dy 


(VI. 11) 


We  now  show  that  eqs . 
Hamilton’ s  canonical  equations 


(5*6)  follow  from  the 


(VI .12) 
(VI. 13) 


siecto  - 


l:  o  :  :  u  '' 
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where 

H  (t)  =  f  3 £  (x,  t)  dx 


(VI .14) 


is  the  total  Hamiltonian  of  the  system,  and  "slash”  denotes 
functional  derivative.  The  functional  dependence  of  is^ 

when  varying  H  with  respect  to  tt. 


It  =  X  (->  If  ) 


so  that  the  usual  relation 


$H  _  d  X  d  ,  d  X 

pi t  ~  d  i t  3x  ^(^7r) 

'-5x' 


(VI. 15) 


must  hold.  With  an  argument  similar  to  that  used  In  deriving 
eq.  (6),  we  can  show  that  the  relation 


d  #  ZX 
^  '  T7W) 

v  ox; 


(vi.i6) 


also  holds. 


The  first  of  the  Hamilton’s  canonical  equations, 
eq.  (12),  now  reads 


bz  a  bX 

Tnr  b(^L) 

loxy 

=  -Jk  /  K(x,y)^(y)  dy  -  •  (-||  gf) 

which  is  easily  seen  to  be  identical  with  eq.  (5). 
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Likewise,  the  remaining  canonical  equation  gives 
rise  to  the  conjugate  Schrftdinger  equation. 


VI. 3  We  are  now  in  a  position  to  apply  the  technique  of 
second  quantization.  Following  the  usual  rule,  we  impose 
the  equal -time  commutation  relation 


[f(x)s  tt(x')]  =  if}  6(x-xt) 

t  t 

Since  it  =  if}if/1  ,  where  ip  1  denotes  the  Hermitian  conjugate  of 
the  operator  if/ ,  we  obtain 


[ip(x)s  *t(x')l  _  8(x-xT)  (VI. 17) 

The  equation  of  motion  of  an  operator  A  with  no  explicit 
dependence  on  time  is  taken  to  be 


ifiA  =  [  A,  H ] 


(VI .18) 


In  particular. 


rh if/  =  [if/,  H] 


Making  use  of  eqs.  (11, l4),  we  find 

00  v,  2 


ifif  =  UJ  Tij  Mr  .  dx'  ]  + 

-00 


(vi.  19) 


[ip,  ff  K(x ' ,  x" )  V  (x '  )  ip(x")  dx’dx"] 
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To  evaluate  the  first  commutator,  we  simplify  the  integral : 


f  dx 


-00 


2m  dx  T  dx  T 


t  _ 


^2,/,+  d^ 

^x7- 


<*>  2  .  ^2, 

11  /  f  t  U_  dx 


-00 


2m 


dx 


2 


tl  r  *f  &2  » 

2m  J  *  9x’2 


dx 


Hence  the  first  commutator  becomes 


i*.  «*■  i 


^  It  ,  I  *f  U  ax'  ] 


2m 


dx 


2 


1h  /  U  (x),H(x’)J  7-^  dx 


dx 


di2  r  R  ,  .  x  d2^  , 

2H  /  6(x-x’)  ^  dx 


dx~ 


1c2  #  d2^ 

2m  %  2 
ox 


The  second  commutator  is 


ff  K(x',x")  [ip  {x)ip^  (x  ’  )ip  (x"  )  -  ^(x'  )^(x")|/(x)  }  dxTdx" 

=  ff  K(x',x")  {^(x  T  )ip(x)ip(x"  )  +  6  (x-x  ’  )ip  (xM  )  -ip^  (x  '  )f{x"  )ip(x)  } 

dx ' dx" 
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=  /  K(x,x")  ^(x")  dx" 
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Hence 


2  -n2 


=  "  +  /  K(x;x’)^  (*')  dx1  , 

chc 


which  is  formally  identical  with  eq.  (5) 


The  Hermit! an  operator 


N  =  /  dx 


(VI. 19) 


is  commonly  assumed  to  represent  the  number  of  particles  in 
the  field.  Its  time  derivative  is  given  by 

iliN  =  [N,H] 

=  [ JVtydx,  /  dx’  ]  +  [  JipTip  dx,  . 


/ /  K(x.T  ,x")^ (x1  )^(x”)dx,dx"  ] 


The  second  commutator  is 


///  dx  dx T  dxM  \  if/^  (x)'if/  (x)^(x '  )f  (x"  )  -i/j^  (x  ’ )  ^  (x"  )^(x)  tfj  (x) 


1 


‘ K (x  T , x" ) 


///dx  dx  T  dx"£^(x)  (x)+6  (x-x (x"  )• 

^  (x  1  )  fif/T  (x)if/  (x"  )  +  6  (x-x"  ^  f  (x)J  *K(x  t  ,  x" 
///  dx  dx  T  dx"  j^^(x)^  (x"  )  5  (x-x  T  )  -^(x  ’  )  ip  (x)  6  (x-xM  )J  K(xfx") 

///  dxT  dx"  ^(xT)^(x,T )  K(x’,x")  -  /  dxT  dx"  /(x')^(x")  K(xT,x") 


0 


iiM 
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(ex. iv) 
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A  similar  calculation  shows  that  the  first  commu¬ 
tator  also  vanishes.  Hence 

N  =  0  ;  (VI. 20) 

therefore  N  is  a  constant  of  motion. 

We  see  that  the  essential  relations  appearing  in 
a  quantized  Schrftdinger  wave  field  for  local  potentials 
can  be  generalized  to  include  nonlocal  interactions. 
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